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THE  COMPRESSION  OF  WET  SNOW 


S.C.  Colbeck,  K.A.  Shaw  and  G.  Lemieux 


INTRODUCTION 

A great  deal  of  information  has  been  generated 
about  the  deformation  of  snow  at  subfreezing  tem- 
peratures (e.g.,  see  review  by  Mellor  (1975)1,  but 
relatively  little  is  known  about  the  deformation  of 
snow  containing  liquid  water.  The  Japanese  have  pro- 
vided some  valuable  information  on  the  compression 
of  wet  snow  (e.g.,  Kinosita  1963,  Ito  1969,  Tusima 
1973,  Wakahama  et  al.  1974,  Wakahama  1975)  show- 
ing that,  with  increasing  liquid  water  content,  snow  is 
more  easily  compressed.  Yarkin  et  al.  (1972)  reported 
that  snow  could  be  strengthened  by  adding  small 
quantities  of  salt,  a fact  that  has  been  used  by  ski  area 
operators  for  hardening  melting  snow  surfaces.  The 
potential  for  the  application  of  knowledge  of  these 
and  other  features  of  the  peculiar  properties  of  wet 
snow  to  improve  snow  removal  from  highways,  control 
wet  snow  avalanches,  and  achieve  greater  mobility  over 
wet  snow  covers  requires  a basic  understanding  of  the 
material  properties  of  wet  snow. 

Many  studies  have  been  made  of  phase  equilibrium 
of  two-  or  three-component  (water,  air  and/or  ionic 
impurities)  systems.  For  exampie.  Berry  (1959)  in- 
vestigated the  dependence  of  the  triple  point  of  water 
on  the  size  of  ice  particles.  Dorsey  (1940)  reviewed 
the  effects  of  pressure,  dissolved  air,  and  dissolved 
ionic  impurities  on  the  liquid-solid  phase  equilibrium 
temperature.  Defay  and  Prigogine  (1951)  placed  the 
thermodynamic  theory  of  multi-phase  and  multi- 
component  systems  in  a convenient  form  for  applica- 
tion to  porous  materials  just  as  Dufour  and  Defay 
(1963)  did  for  applications  to  clouds.  Using  similar 
methods,  Colbeck  (1973)  developed  the  relationships 
describing  phase  equilibrium  in  wet  snow,  showing  how 
the  liquid  water  content  and  grain  size  affected  the 
delicate  thermodynamics  of  the  material. 

Wet  snow  is  a very  unusual  material  in  that  its 
response  to  an  applied  load  can  be  dominated  by  small 
temperature  differences  between  the  stressed  and 


stress-free  surfaces  of  the  individual  particles.  These 
temperature  differences  arise  because  the  equilibrium 
temperature  at  a solid/liquid  interface  decreases  with 
increasing  pressure;  hence,  the  stressed  contact  surfaces 
between  neighboring  particles  exist  at  a temperature 
which  is  generally  lower  than  the  particles’  stress-free 
pore  surfaces.  Colbeck  (1976)  showed  that  these  tem- 
perature differences  around  the  surfaces  of  individual 
particles  account  for  the  rapid  deformation  of  wet 
snow,  including  an  explanation  of  the  effect  of  liquid 
water  content  on  the  compressibility  of  the  snow. 
Colbeck  and  Parssinen  (in  press)  expanded  this  model 
of  compressibility  to  include  dissolved  ionic  impurities 
and  dissolved  air  but  calculated  rates  of  deformation 
which  were  generally  too  large. 

In  the  present  study,  simultaneous  grain  growth, 
the  double  charged  layer  at  phase  boundaries,  and  an 
additional  mechanism,  power  law  creep,  are  included 
in  the  model.  Also,  a long  series  of  experiments  have 
been  completed  and  are  reported  here  along  with  a 
complete  development  of  the  model  of  the  compressi- 
bility of  wet  snow. 

EXPERIMENTAL  PROCEDURES 

Samples  of  wet  snow  were  prepared  by  sieving  old 
dry  snow  and  immersing  the  sieved  snow  in  liquid 
water  for  about  24  hours  before  loading  the  samples 
into  the  experimental  apparatus.  The  mixture  was  then 
packed  in  a test  cylinder  (see  Fig.  1 ),  using  a careful 
procedure  in  order  to  achieve  a uniform  packing  of  the 
samples  with  a density  of  about  0.48  and  an  average 
grain  diameter  of  about  1 mm.  These  samples  were 
loaded  with  a “dead  weight"  whose  movement  was 
measured  with  a displacement  transducer  (LVDT)  in 
order  to  calculate  the  ice  density  (ice  mass  per  unit 
volume)  throughout  the  experiment.  The  most  diffi- 
cult part  of  the  experimental  procedure  was  controlling 
the  temperature  of  the  sample  by  maintaining  the 
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Figure  I . A diagram  of  the  compression  apparatus 
located  inside  a temperature  controlled  cabinet. 

ice-water  bath  at  the  same  temperature  as  the  pore 
water  in  the  sample  during  the  two-week  duration  of 
each  experiment.  In  order  to  achieve  proper  tempera- 
ture control,  the  ice-water  bath  had  to  be  maintained 
at  the  same  salinity  as  the  water  in  the  sample  and 
had  to  be  maintained  with  approximately  equal  parts 
of  ice  and  water.  This  required  frequent  testing  of  the 
sample  and  bath  salinities  and  additions  of  salt  and/or 
ice  to  the  bath.  Many  experiments  were  rejected  be- 
cause of  evidence  of  melting  in  the  sample  or  freezing 
around  the  weight  upon  completion  of  an  otherwise 
successful  run. 

Two  sets  of  experiments  were  successfully  com- 
pleted showing  the  effect  of  changing  the  sodium 
chloride  content  and  the  capillary  pressure  (results  are 
shown  later).  The  capillary  pressure  was  held  constant 
during  each  experiment  by  setting  the  bottom  of  the 
water  filled  tube  at  a fixed  position  below  the  bottom 
of  the  sample.  Assuming  the  "liquid  path”  is  continu- 
ous, this  procedure  sets  the  "tension"  in  the  liquid 
water  in  the  sample  at  a value  which  is  constant  through- 
out the  experiment.  As  the  compression  proceeds,  the 


sample  discharges  water  to  maintain  the  imposed 
capillary  pressure.  Unlike  other  experiments  reported 
in  the  literature,  there  was  no  "saturated  layer”  at  the 
bottom  of  these  samples  and,  although  the  liquid  water 
content  varied  in  all  experiments,  the  parameter  of 
thermodynamic  significance  (capillary  pressure)  re- 
mained constant  during  each  run. 

GRAIN  GROWTH 

Although  the  snow  grains  were  immersed  in  liquid 
water  for  24  hours  before  being  packed  into  the  com- 
pression device,  grain  growth  continued  throughout 
the  duration  of  the  experiments.  After  being  immersed 
for  24  hours,  the  snow  grains  achieved  an  average  diam- 
eter of  about  1 mm,  and  after  being  stressed  and  im- 
mersed for  2 weeks  they  reached  an  average  diameter 
of  about  2 mm.  In  similar  experiments,  Wakahama 
(1968)  measured  the  distribution  of  sizes  of  snow 
grains  shown  in  Figure  2.  From  Figure  2,  we  calculate 
the  average  radius  z using 


and  the  ratio  of  particle  numbers  using  conservation  of 
mass,  or 

Nr- 3 = yVoz3  . (2) 

These  results  are  given  in  Table  I for  the  distributions 
shown  in  Figure  2. 

Working  with  bismuth,  which  has  the  typical  melting 
temperature-radius  relationship,  Peppiatt  (1975)  showed 
that  the  number  of  bismuth  particles  in  the  presence  of 
the  melt  decreases  exponentially  as  a function  of  time 
and  chose  to  represent  his  observations  by  the  relation- 
ship 

N = Nq  e-^h  . (3) 

From  eq  2 it  follows  that  we  could  represent  the 
average  particle  radius  by 

z = Zq  . (4) 

However,  this  suggests  that  the  rate  of  increase  of 
particle  radius  increases  with  time.  Wakahama’s  (1968) 
data  and  our  data  clearly  show  that  the  rate  of  increase 
of  z decreases  with  time;  hence,  we  adopt  the  expression 

z = Zo  + 1.582  (z,-Zo)(1 -e-f/^)  (5) 
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Gram  Diameter,  mm 

Figure  2.  The  distribution  of  sizes  of  snow  grains  initiaiiy,  immersed  in  water  for  three  days,  and 
immersed  in  water  for  six  days  (after  Wakahama  1968).  We  smoothed  the  upper  end  of  the  dis- 
tribution for  six  days. 


I.582(r,.-r,)(l-e' 


(o)  Wokohama  (toturated,  no  load) 
(^)  Colbeck  (soturotod,  looded) 


0 0.5  I0«I0* 

Time  (s) 

Figure  3.  The  increase  in  grain  radius  shown  with  time  for  r^  = 1 mm  and 
T = 1,209,600  s. 


Table  I.  Calculations  from  distributions  of  Wakahama  (1968). 


Figure  4.  The  relative  increase  in  grain  radius  and  decrease  in  number  of  particles 
shown  against  the  characteristic  time  for  r^  = 2rQ. 


which  is  shown  in  Figure  3 together  with  the  experi- 
mental data.  It  follows  that 

N = /Vo[1  -r  1.582(r,-ro)(1  -e-'/n/rol'^  (6) 

and  that  drjdt  decreases  exponentially  with  time  (see 
Fig.  4)  as  shown  by  all  of  the  existing  experimental 
data. 

There  are  two  reasons  why  grain  growth  must  be 
considered  simultaneously  with  compression.  First, 
both  processes  are  heat  dominated  and  the  heat  source, 
refreezing  on  the  stress-free  surfaces  of  the  larger 
particles,  must  be  shared  between  the  stressed  surfaces 
and  the  shrinking  particles.  Second,  the  fact  that  seven 
out  of  eight  of  the  particles  melt  away  during  the  two 
weeks  of  compression  indicates  that  many  vacancies 
must  temporarily  exist  in  the  packing  as  particles  melt 
away  from  their  neighbors.  Both  Wakahama’s  (1 975) 
observations  and  Colbeck’s  (1973)  calculations  show 
that  once  a particle  begins  melting  it  disappears  quickly. 
Thus  only  a small  fraction  of  the  particles  is  shrinking 
at  any  time  but,  as  shown  later,  these  particles  reduce 
the  average  number  of  contacts  per  particle  and  in- 
crease the  average  stress  per  contact  accordingly. 


rather  than  a systematic  packing  (the  effect  of  this 
number  on  the  rate  of  compression  is  shown  later). 
Because  the  shrinking  particles  are  under  negligible 
stress  at  their  contacts,  the  average  number  of  stressed 
contacts  per  particle  n is  given  by 

n = 6-^2NJN  (7) 


where  NJN  is  the  fraction  of  the  particle  population 
which  is  shrinking  and  the  coefficient  (12)  arises  be- 
cause one  shrinking  particle  removes  12  contacts  (6  of 
its  own  and  6 from  neighboring  particles).  The  rate  of 
change  of  the  numbei  of  particles  is  given  by 


dNIdt  = -NJtQ. 
Flence 


dt 


or 


(8) 


(9) 


n = e-51-^ ^^2 . (10) 

Ted-'  /-o  + 1.582(/-,-/-o)(1 -e-'/’-) 


STRESSED  PARTICLE  CONTACTS 

Visscher  and  Bolsterli  (1972)  found  an  average  of 
6.4  contacts  per  sphere  in  random  packings  of  equal 
size  particles.  For  simple  cubic  packings  there  are  6 
contacts  per  particle,  which  we  use  here  as  the  upper 
limit  of  the  average,  although  wet  snow  is  a random 


A typical  case  (see  the  parameters  listed  in  Table  II) 
where  n increases  with  time  is  shown  in  Figure  5.  Over 
10®  %,n  increases  from  about  2 to  5,  while  the  stress 
at  the  bonds  decreases  rapidly.  In  this  example  about 
two-thirds  of  the  particles  are  stressed  initially;  hence  the 
average  number  of  stressed  bonds  per  stressed  particle 
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Figure  5.  Computed  values  of  the  average  number  of  contacts  per  particle  and  the 
fluid  stress  at  a bond  shown  against  time  for  parameters  shown  in  Table  II. 


Table  II.  Parameters  for  various  figures. 


Figure 

"bulk 

(Pa) 

<■0 

(mm) 

'r 

(mm) 

r 

(s) 

'd 

(s) 

Pc 

5 

11,790 

0.0007 

0.5 

1 

1 ,209,000 

86,400 

0 

8 

11,790 

0.0007 

0.5 

1 

1,209,000 

86,400 

0 

9 

11,790 

0.0007 

0.5 

1 

1,209,000 

86,400 

0 

10 

10® 

0.0007 

0.5 

1 

1,209,000 

86,400 

0 

11 

4,790 

0.0007 

0.5 

1 

1 ,209,000 

86,400 

4 values 

12 

11,790 

3 values 

0.5 

1 

1,209,000 

86,400 

0 

13 

11,790 

4 values 

0.5 

1 

1,209,000 

86,400 

0 

14 

4,790 

0.0007 

3 values 

4 values 

1 ,209,000 

86,400 

0 

15 

4,790 

0.0007 

0.5 

1 

1,209,000 

3 values 

0 

16 

4,790 

0.0007 

0.5 

3 values 

1,209,000 

86,400 

8,923  Pa 

17 

4,790 

0.0007 

0.5 

1 

1,209,000 

86,400 

0 

18 

4,790 

0.0007 

0.5 

1 

1,209,000 

86,400 

0 

19 

4,790 

' varies 

2 values 

2 values 

1 ,209,000 

86,400 

varies 

is  about  3 initially.  Thus  there  are  enough  stressed 
contacts  to  make  the  particles  stable  geometrically  at 
all  times. 


FLUID  PRESSURE 

As  the  densification  proceeds,  the  thickness  of  the 
liquid  film  between  the  particles  approaches  10*^  mm; 
hence  we  must  consider  the  electrostatic  forces  which 
arise  due  to  the  double  charged  layers  at  each  solid-liquid 
boundary  (see  cover  photo).  It  is  generally  agreed  that 
these  forces  decay  exponentially  with  separation 


(Israelachvili  and  Ninham  1977,  Aveyard  and  Haydon 
1973,  p.  43,  Adamson  1967,  p.  218)  and  increase  with  de- 
creasing salinity  (Israelachvili  and  Adams  1976).  While 
the  van  der  Waal  forces  at  this  separation  are  small,  the 
repulsive  forces  between  the  double  layers  are  signifi- 
cant and  counteract  the  applied  load,  especially  at  low 
salinities.  From  the  experimental  results  of  Israelachvili 
and  Adams  (1976),  who  used  KNO3  solution,  we  de- 
duce that  the  repulsive  force  per  unit  area  between  two 
flat  surfaces  separated  by  a NaCI  solution  is 

exp  (-2.3  X 1 0®  yjs^  d)  (units  of  Pascals) 
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where  d is  in  units  of  millimeters  and  5),  in  molality. 
The  coefficient  (0.42)  is  about  one-fourth  that  of 
Israelachvili  and  Adams  because  NaCI  is  a stronger 
structure  former  than  KNO3  (Frank  and  Tsao  1954). 

The  average  stress  at  an  interparticle  contact  de- 
creases rapidly  after  the  stress  is  applied  because  of 
the  increasing  number  of  contacts  available  to  bear  the 
load  and  the  increasing  size  of  the  contacts  as  pressure 
melting  occurs.  The  increasing  number  of  contacts 
decreases  stress  by  a factor  of  approximately  6ln, 
while  the  increasing  area  decreases  the  stress  by  the 
ratio  of  the  area  of  the  particle  to  the  area  of  the  bond. 
Accordingly,  the  stress  in  the  liquid  at  the  bond  is 
given  by 


PACKING  GEOMETRY 

The  depth  of  the  melted  cap  D and  the  radius  of 
the  bond  can  be  calculated  at  each  stage  of  the  de- 
formation from  the  conservation  of  the  total  ice  mass 
in  the  sample,  or  (see  Fig.  6) 

Srg/Vj,  = 8r3/V-n/VA(A2  + 3/-g)  (12) 


where 

A = /■p-(/-2-r2)14  . (13) 

Given  these  values  we  can  calculate  the  ice  density  of 
the  sample  from 

which  is  an  exact  expression  for  simple  cubic  packings 
where  n is  an  integer.  The  computations  are  made 
iteratively  by  calculating  the  rate  of  increase  of  D at 
each  step  and  then  updating  the  geometrical  variables. 


TEMPERATURE  DISTRIBUTION 
AND  HEAT  FLOW 

The  rate  of  increase  of  D is  given  by  the  heat  flow 
from  the  stress-free  to  the  stressed  surfaces  minus  a 
correction  for  the  heat  diverted  to  the  shrinking 
particles.  Accordingly, 


Figure  6.  The  cross-sectionai  view  of  a typicai  pvticie 
symmetricaiiy  deveioping  four  melt  caps.  The  smaller 
circle  represents  the  Initial  particle  radius  and  the  larger 
circle  represents  the  Increased  particle  radius  due  to 
the  refreezIng  on  the  stress-free  surfaces. 


dt  Lp\l  2ti  ^ 


(15) 


where 


2ri-6r^h 


H 


2ri-6r„h 


(16) 


L£_rA:  if  “S  '~P"  < 1 

H H 


and  AT  is  the  temperature  difference  between  the 
stressed  and  stress-free  surfaces.  The  first  term  in  eq 
1 5 represents  the  heat  flow  due  to  the  temperature 
difference  AT,  whereas  the  second  term  reduces  the 
heat  flow  to  account  for  the  fact  that  simultaneous 
particle  growth  diverts  some  heat  away  from  the  only 
source,  refreezing  on  the  stress-free  surfaces  of  the  ice 
particles.  The  first  term  takes  the  varying  number  of 
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stress  contacts  per  particle  n into  account  because  the 
average  path  length  for  heat  flow  decreases  with  in- 
creasing/?. Also,  the  ratio  of  the  surface  areas  acting 
as  heat  sources  and  heat  sinks  is  accounted  for  by  A. 
The  second  term  is  derived  from  a total  heat  balance. 

The  temperature  difference  AT  is  equal  to  the  phase 
equilibrium  temperature  of  the  stress-free  surface 
minus  the  phase  equilibrium  temperature  of  the  stressed 
surface  minus  a small  correction  for  supercooling  on 
the  stress-free  surface, 


There  are  two  regimes  of  liquid  saturation,  pe/?dt//<7r 
and  funicular,  which  must  be  considered  separately. 

In  the  pendular  regime  of  lower  saturation  where  the 
air  phases  exist  in  continuous  paths  throughout  the 
pore  space,  the  phase  equilibrium  at  the  stress-free 
surfaces  is  (Colbeck  1973): 


Tq  ^Tq  ^ 

'■p 


(18) 


where  the  first  term  arises  due  to  capillary  pressure, 
the  second  due  to  the  size  of  the  ice  particles,  and  the 
third  due  to  dissolved,  ionic  impurities.  At  higher 
liquid  saturations  where  the  liquid  phase  exists  in  con- 
tinuous paths  throughout  the  pore  space,  the  phase 
equilibrium  temperature  at  the  stresvfree  surfaces  is 


■ - ^^0  /_L  J_\££i  iZo.  ^^0  /c 
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(19) 
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Figure  7.  Two  particles  pushed  to- 
gether melt  at  the  stressed  surfaces. 

The  melt  water  Is  removed  in  a thin 
film  and  re  freezing  takes  place  on 
the  stress- free,  rounded  surfaces. 

temperature  decrease  due  to  the  presence  of  dissolved 
ionic  impurities.  Note  that  no  pressure  melting  would 
occur  if  water  were  an  ordinary  substance  and  con- 
tracted upon  freezing  because  the  first  term  would  then 
be  positive  and  the  solid-liquid  equilibrium  temperature 
would  increase  with  increasing  stress. 

FLUID  AND  IMPURITY  FLOWS 


where  the  first  term  arises  due  to  the  radius  of  the  The  concentrations  of  dissolved  air  and  ionic 

trapped  air  bubbles.  The  temperature  of  the  stressed  impurities  5|,  in  the  melt  water  being  flushed  through 

surfaces  at  the  inter-particle  conticts  (see  Fig.  7)  is  the  thin  film  between  the  particles  is  of  critical  impor- 

given  by  tance  to  the  rate  of  compression  of  wet  snow.  Drake 

and  Shreve  (1973)  showed  that  wires  moving  through 
r - ^0/1  1 \ ice  by  regelation  behave  similarly,  although  the  source 

•>  ~ T~  ~ ^ ^ \ C,  / ~L  of  the  dissolved  impurities,  diffusion  from  the  pore 

water  into  the  very  thin  film  of  liquid  between  the 
particles  (see  Fig.  8),  would  seem  to  be  very  slow  be- 
(20)  cause  of  its  small  size.  Despite  this  small  size,  the 
particle  path  for  diffusion  is  also  very  small  and  the 
where  the  first  term  is  the  temperature  decrease  due  to  concentration  gradients  are  fairly  large, 

the  stress  exerted  on  the  solid  by  the  liquid  (the  double  The  melt  water  moves  radially  outward  through  the 

layer  forces  cancel  at  each  interface;  hence  Ty,  is  deter-  liquid  film  at  a rate  given  by  (Moore  1965) 

mined  by  fluid  pressure  only),  the  second  term  is  the 


0.4  0.6 

Time  (s) 


Figure  8.  The  computed  size  of  the  liquid  film  between  the  particles  and  the  fraction 
of  dissolved  impurities  reaching  the  bond  shown  against  time  for  typical  conditions. 


Continuity  requires 

V = 

d Pj  dt 

Hence  the  film  thickness  is 


Here  we  assume  distinct  phase  boundaries  and  a film 
of  homogeneous  liquid  whose  flow  is  not  affected  by 
charge  separation.  Some  modifications  to  these  as- 
sumptions are  suggested  later. 

Assuming  a quasi-steady  balance  between  the  move- 
ment of  the  impurities  with  the  melt  water  and  the 
diffusion  of  the  impurities  against  the  flow  of  the 
melt  water, 


vC*D,^  = 0 
' dr 


or  the  distribution  of  the  impurities  is  given  by 


where  C),  is  the  concentration  at  the  center  of  the 
contact.  Assuming  that  the  contact  surfaces  arc  flat 


[d  * f(r)l , we  can  simplify  eq  25,  but  there  is  a 
complicated  feedback  involving  the  impurity  concen- 
tration, the  double  charged  layer,  the  temperature 
difference,  the  rate  of  melting,  the  melt  water  dis- 
charge rate  and,  again,  the  impurity  concentration. 

This  feedback  loop  demonstrates  the  delicate  nature 
of  the  pressure  melting  mechanism. 

Another  mechanism,  nonlinear  viscous  creep  of  the 
ice  in  the  region  of  the  stressed  contacts,  is  also  con- 
sidered. Wilkinson  and  Ashby  (1975)  show  that  for  a 
material  with  a power  law  relationship  between  stress 
and  strain  rate,  or 

e = Ao^  (26) 

the  rate  of  densification  (dp/dt)  due  to  the  creep  of  the 
material  at  the  contact  is 

dt  4 fr  ^2  _2  m / • 


Here  A f is  4 r^  for  a simple  cubic  packing.  We  take  m 
as  three  and,  from  the  results  of  Colbeck  and  Evans 
(1973),  we  take  A equal  to  1.28x10*23  (Pa-m  5-I)  for 
polycrystalline  ice  in  thermal  equilibrium  with  the  melt. 
For  a typical  case,  calculated  values  of  the  rate  of 
densification  are  shown  for  both  the  total  and  non- 
linear viscous  deformation  mechanisms  in  Figure  9. 
Clearly  regelation  is  the  dominant  mechanism  in  this 
example,  but  because  the  rate  of  densification  due  to 
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Figure  9.  The  computed  rates  of  densification  shown  for  the  total  and  the  nonlinear 
viscous  deformations. 


power  law  creep  increases  as  a high  power  of  stress, 
whereas  the  temperature  depression  causing  regelation 
increases  as  only  the  first  power  of  stress,  the  nonlinear 
mechanism  is  relatively  more  important  at  higher 
stresses. 

The  ice  density  is  shown  as  a function  of  time  in 
Figure  10  for  an  applied  stress  of  10^  Pa.  This  shows 
that  even  at  stresses  as  large  as  the  overburden  pressure 
at  the  firn-ice  transition  in  a temperate  glacier,  regela- 
tion is  still  the  dominant  mechanism.  This  conclusion 
is  reinforced  by  observations  le.g.,  Vallon  et  al.  (1976)) 
which  have  shown  a correlation  between  the  rate  of 
firniflcation  and  the  liquid  water  content  of  the  snow. 


This  correlation  can  .k  readily  explained  in  terms  of 
the  regelation  described  here. 


RESULTS 


As  already  shown  by  the  results  given  in  Figures  5, 

8, 9 and  10,  a model  has  been  constructed  using  the 
concepts  described  in  the  previous  section.  The  cxper- 
mental  results  at  various  liquid  water  contents  are  shown 
together  with  computed  curves  in  Figures  1 la  and  11b. 
Both  sets  of  curves  clearly  show  the  decreased  rate  of 
densification  at  lower  liquid  contents  (higher  capillary 
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Figure  10.  The  ice  density  shown  against  time  for  regelation  plus  nonlinear 
viscous  creep  and  regelation  alone.  The  difference  between  the  two  lines  is 
the  contribution  of  the  power  law  creep. 


pressures)  resulting  from  a reduced  temperature  dif- 
ference as  explained  by  eq  18.  There  are  consistent 
differences  between  the  experimental  and  computed 
results  which  probably  arise  from  the  approximations 
made  about  heat  flow  and  packing  geometry,  some 
uncertainty  over  the  exact  value  of  the  initial  density 
in  the  individual  experiments,  and  some  uncertainty 
over  the  continuity  of  the  “liquid  path"  in  the  sample 
with  the  largest  value  of  capillary  pressure. 

Similarly,  Figure  12  shows  experimental  and  com- 
puted results  for  snow  saturated  at  various  salinities. 

The  decreased  rate  of  densification  with  increasing 
salinity  is  apparent  in  both  sets  of  curves,  although 
there  are  differences  between  the  computed  and 
measured  results.  Given  the  severe  experimental  limi- 
tations inherent  in  long-term  loading  of  samples  of 
wet  snow,  we  do  not  know  whether  the  major  prob- 
lems are  with  the  fine  points  of  the  experiments  or 
the  model.  Because  of  the  severe  difficulties  of  working 
with  this  material  and  the  variety  of  parameters  which 
are  involved,  the  model  is  very  useful  to  help  under- 
stand the  physical  processes,  to  identify  important 
parameters,  and  to  make  quick  estimates  of  the  com- 
pressibility under  widely  varying  conditions. 

Figure  1 3a  shows  the  wide  range  of  responses  com- 
puted for  different  salinities.  The  largest  change  occurs 
as  salinity  decreases  to  10'^  where  the  double  charged 
layer  begins  to  dominate  other  processes.  At  lower 
salinities  the  double  charged  layer  decreases  the  stress 


in  the  thin  liquid  film,  thus  reducing  the  temperature 
gradient  and  heat  flow  to  the  melting  surfaces.  This 
effect  of  the  double  layer  is  clearly  demonstrated  by 
comparing  Figures  13a  and  13b.  In  Figure  13b  where 
the  electrostatic  forces  were  ignored,  the  computed 
densities  at  low  salinities  are  much  higher  than  those 
observed  experimentally.  Thus,  the  importance  of  the 
double  layer  at  low  salinities  is  established. 

Sodium  chloride  has  a small  effect  on  the  rate  of 
growth  of  snow  grains,  but  this  effect  is  small  com- 
pared with  the  greatly  reduced  rate  of  growth  at  low 
liquid  water  contents  (Wakahama  1975).  No  doubt 
the  imposed  load  also  changes  the  rate  of  grain  growth; 
hence  it  is  necessary  to  investigate  these  matters  here. 
Figure  14  shows  the  computed  increase  in  ice  density 
with  time  for  four  situations  with  different  grain  radii. 
For  a constant  grain  radius  of  0.5  mm,  the  largest  rate 
of  densification  occurs  because  there  is  no  heat  diversion 
to  growing  particles  and  because  the  heat  flow  path 
length  is  smallest  for  smaller  particles.  At  a constant 
radius  of  1 mm,  the  computed  rate  of  densification  is 
reduced  markedly  because  of  the  larger  path  length  for 
heat  flow.  The  measured  increase  in  ice  density  with 
time  for  a particle  size  of  3 to  4 mm  shown  in  Figure 
14  confirms  this  idea  about  the  heat  flow  path  length. 

It  is  rather  suprising  that  there  is  little  difference  be- 
tween the  densification  of  particles  of  a constant  1-mm 
radius  and  those  which  grow  from  0.5  to  1 mm  during 
the  densification.  The  densification  of  particles  which 
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b.  Computed  results. 

Figure  / /.  ke  density  shown  against  time  for  various  capillary  pressures  for  both 
experimental  (a)  and  computed  (b)  results 
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Figure  12.  Ice  density  shown  against  time  for  various  sodium  chloride  contents 
(g-moles  Nad  per  kilogram  of  solution  times  dissociated  ions  per  moiecule)  for 
experimental  and  computed  results. 
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Figure  13a.  Computed  ice  density  shown  against  time  for  various  salinities 
(g-moies  per  kilogram  of  solution  times  dissociated  ions  per  molecule). 
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grow  from  0.5  to  0.75  mm  during  the  densification  is 
intermediate  between  particles  which  remain  at  0.5 
mm  and  those  which  grow  to  1 .0  mm.  Clearly,  the  rate 
of  particle  growth  has  a major  effect  on  the  rate  of 
densification  in  saturated  snow  and  the  characteristics 
of  the  growth  must  be  well  known  in  order  to  quantify 
the  compression  of  wet  snow. 


The  initial  particle  size  is  imposed  by  the  situation, 
whereas  the  final  particle  size  is  partly  the  result  of  the 
experimental  parameters.  All  of  the  required  informa- 
tion can  be  supplied  from  knowledge  of  the  conditions 
under  which  the  deformation  takes  place  and  observa- 
tions of  particle  growth  under  similar  conditions.  The 
average  life  expectancy  of  a shrinking  particle  Zq  can 
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Figure  14.  Four  examples  of  computed  ice  density  are  shown  against  time  for 
various  particle  radii  and  rates  of  particle  growth.  One  experimental  result  Is 
shown  for  large  ice  chips  of  3 to  4-mm  diameter. 


Figure  IS.  Computed  ice  density  shown  against  time  for  various  life  expectancies 
for  an  average  drinking  particle. 


be  determined  using  techniques  like  Wakahama’s 
(1975)  time-lapse  photography  or  Colbeck’s  (1973) 
calculations  of  particle  capture.  Here  we  have  chosen 
a value  of  1 day  (86,400  s)  which  is  in  the  range  of  the 
life  expectancies  observed  and  calculated.  As  shown  in 
Figure  1 5,  the  computed  rate  of  densification  is  not  too 
sensitive  to  this  parameter;  hence  it  is  not  necessary  to 
adjust  the  life  expectancy  to  accommodate  changes  in 
the  other  parameters  in  order  to  use  this  model. 

At  large  capillary  pressures  or  large  salinities,  the 
effect  of  simultaneous  particle  growth  is  greatly  re- 
duced because  of  the  altered  temperature  field.  Ap- 
parently, simultaneous  heat  diversion  to  growing 


particles  is  not  of  much  consequence  when  the  tem- 
perature of  the  stress-free  particle  surfaces  is  reduced 
by  low  liquid  contents  or  large  concentrations  of  im- 
purities in  the  pore  water.  Figure  16  shows  the  effect 
of  changing  the  particle  radius  after  growth  r^  for  an 
initial  particle  radius  of  0.5  mm  at  a capillary  pressure 
of  8923  Pa.  It  is  clear  from  Figures  14  and  16  that 
simultaneous  particle  growth  has  a greatly  reduced 
effect  on  densification  at  large  capillary  pressures; 
hence  the  varying  rate  of  particle  growth  with  liquid 
water  content  is  not  the  important  parameter  we  might 
have  supposed.  This  greatly  simplifies  the  use  of  this 
model,  since  for  most  purposes  we  can  ignore  the 
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Figure  16.  Computed  ice  density  shown  against  time  for  various  values  of  final 
particle  radius  at  a large  capillary  pressure.  A final  value  of  0.5  represents  no 
particle  growth,  whereas  a final  value  of  1.0  mm  represents  a doubling  of  the 
average  particle  size. 


Figure  / 7.  Computed  ice  density  shown  against  time  for  two  values  of  the  upper 
limit  of  the  average  number  of  contacts  per  particle. 


variation  of  particle  growth  with  liquid  water  content 
and  salinity. 

We  used  a simple  cubic  packing  for  assistance  in 
conceptualizing  at  several  steps  in  previous  sections; 
for  example,  in  the  derivation  of  eq  7,  we  used  a maxi- 
mum of  6 contacts  per  particle  rather  than  the  average 
of  6.4  calculated  by  Visscher  and  Bolsterli  (1972).  We 
are  now  in  a position  to  show  that  this  assumption  has 
a small  effect  on  the  values  computed  using  this  model. 
Figure  1 7 shows  the  result  of  changing  the  average 
number  of  contacts  per  particle  from  6 to  6.4  and,  as 
with  the  various  values  for  Zq  shown  in  Figure  IS, 
there  is  very  little  error  resulting  from  the  use  of  this 
approximation. 

The  results  shown  in  Figure  9 suggest  that  the  rate 
of  density  increase  is  proportional  to  time  to  the  power 
negative  0.8.  Therefore  we  may  express  the  density  of 
wet  snow  as 


P = ^(®bulk-  Pc<  '■q)  Po  ■ (28) 

This  expression  is  compared  with  the  computed  and 
measured  values  of  density  and  time  for  a particular 
case  in  Figure  18,  where  it  can  be  seen  that,  as  long  as 
the  parameters  can  be  included  properly,  this  simple 
expression  could  be  a highly  accurate  substitute  for  the 
computations.  For  many  purposes,  it  may  be  better  to 
work  with  this  expression,  adjusting  the  coefficient  as 
required  to  account  for  the  particular  parameter  under 
consideration. 

just  as  St.  Lawrence  (unpublished)  showed  for  low 
density  snow  at  subfreezing  temperatures,  it  is  possible 
to  describe  the  deformation  of  high  density,  wet  snow 
by  physical  considerations  on  the  scale  of  individual 
particles.  The  most  significant  difference  between 
subfreezing  and  wet  snow  is  the  regelation  mechanism 
described  here.  For  example,  dry  snow  and  saturated 
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Figure  18.  Computed  and  measured  ice  densities  shown  against  time  and  com- 
pared with  an  aigebraic  expression. 
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Figure  19.  Measured  ice  density  shown  against  time  for  dry  snow  at  -TC 
and  saturated  snow. 


snow  compressed  at  the  same  stress  are  shown  in  Figure 
19.  Upon  application  of  the  load,  the  dry  snow  quickly 
compresses  up  to  a density  of  0.474,  the  approximate 
density  at  which  the  wet  snow  samples  were  packed. 
Above  that  density  the  deformation  of  dry  snow  is 
limited  by  the  slow  processes  of  creep  at  the  stressed 
contacts  and  other  sintering  mechanisms  described  by 
Hobbs  (1965). 

Sintering  in  the  presence  of  the  liquid  phase  in 
metals  is  dominated  by  the  “solution  - deposition” 
mechanism  (Kingery  1958,  Courtney  1977)  where 
smaller  particles  disappear  because  of  their  greater 
solubility.  Like  pressure  melting  in  wet  snow,  solution- 
deposition  is  also  sensitive  to  simultaneous  particle 
growth  (Cannon  and  Lenel  1953).  Also,  sintering  in 
metals  is  enhanced  by  deformation  at  the  stressed 
contacts  (Murray  et  al.  1958,  Wilkinson  and  Ashby  1975). 


THE  LIQUID  FILM 

While  the  concept  of  compaction  by  pressure  melting 
is  well  established  for  wet  snow,  several  limitations  of 
the  model  are  important.  The  simplified  treatment  of 
heat  flow  and  reference  to  a simple  cubic  packing  have 
aiready  been  noted,  but  the  most  serious  limitations  of 
this  model  may  be  the  assumptions  made  about  the 
very  thin  liquid  film  between  the  stressed  particles. 
There  are  several  reasons  for  concern.  First,  Figure  8 
shows  a typical  thickness  of  this  film  of  only  10'*  mm. 
Although  Moore  (1965)  states  that  no  molecular  influ- 
ence can  exist  at  1.27x1  O'*  mm  from  the  surface  in 
films  of  mineral  or  vegetable  oil,  and  Buckley  (as  re- 
ported by  Dorsey  (1940)]  found  no  stationary  film  of 
oil  at  thicknesses  as  large  as  3x10'*  mm,  Derjaguin  |as 
reported  by  Dorsey  (1940))  found  some  elastic 
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"rigidity”  in  water  films  less  than  15x10'^  mm. 
Extrapolating  Derjaguin’s  values  to  10'*,  we  find  that 
viscosity  is  nearly  three  orders  of  magnitude  greater 
than  rigidity.  Nevertheless,  it  is  clear  that  the  effective 
viscosity  of  a film  of  lO'^  mm  thickness  between  two 
ice  surfaces  is  less  than  that  of  homogeneous  water. 
Terzaghi  (as  reported  by  Moore  (1965)|  suggested  an 
equivalent  viscosity, 

p'  = ju(1  + constant/cf)  (29) 

which  would  cause  a rapid  decrease  in  the  discharge 
with  decreasing  film  size.  This  type  of  relationship 
would  be  included  in  the  current  model  if  the  appro- 
priate value  could  be  attached  to  the  constant.  Needs 
[as  reported  by  Moore  (1965)]  showed  that  d has  a 
lower  limit  close  to  molecular  dimensions  when  liquid 
is  displaced  from  between  approaching,  parallel  plates. 
However,  Terzaghi’s  expression  is  adequate  for  our 
purposes  since  the  melting  ice  surfaces,  which  are 
sources  of  the  displaced  liquid,  may  prevent  such  a 
thin  film  from  developing. 

Given  a thin  film  of  water  which  is  partly  hetero- 
geneous, the  rate  controlling  diffusion  of  dissolved  air 
and  ions  through  the  fiim  must  be  considered  carefully. 
Dissolved  nitrogen  molecules  have  a diameter  of  3.1 5 x 
10"^  mm  or  over  3%  of  the  minimum  thickness  of  the 
liquid  film.  Dissociated  ions  Na'*’  and  Cl”  have  diame- 
ters of  0.97x10'^  mm  and  1.81  xIO'^  mm  respectively. 
Given  the  size  of  the  fiim  and  the  fact  that  much  of  it 
may  be  heterogeneous  due  to  the  presence  of  groups 
of  molecules,  it  follows  that  diffusion  may  occur  more 
slowly  in  these  thin  films  than  in  homogeneous  water. 

If  so,  the  diffusion  of  impurities  through  the  film 
would  be  less  than  calculated,  and  the  temperature 
gradient  and  heat  flow  would  be  reduced,  thus  reducing 
the  rate  of  melting  and  the  liquid  flow  through  the  film. 
However  reducing  the  flow  would  increase  the  diffusion 
of  impurities;  hence  both  the  increased  viscosity  and 
reduced  diffusivity  have  limited  effects  because  of  this 
stable  feedback  mechanism  in  the  regelation  process. 

Preferential  incorporation  of  cations  or  anions  with 
the  freezing  liquid  on  the  stress-free  surfaces  can  cause 
charge  separation  in  the  pore  water  (Cross  1972).  This 
may  or  may  no^  “nhance  the  perferential  diffusion  of 
ions  into  the  thin  ^Im,  but  the  double  layer  already 
described  is  clearly  important.  The  treatment  of  that 
double  layer  could  be  improved  if  there  were  direct 
evidence  of  its  strength  for  an  ice  surface  and  a sodium 
chloride  solution. 

CONCLUSIONS 

Wet  snow  is  a very  unusual  material  which,  at  inter- 
mediate to  high  densities,  deforms  predominantly  by 
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regelation.  The  physical  model  of  this  process  re- 
ported here  allows  us  to  investigate  the  effects  of  grain 
growth,  liquid  saturation  and  dissolved  impurities.  It 
is  well  known  that  the  compressibility  of  wet  snow  de- 
creases with  increasing  salinity  and  decreasing  liquid 
water  content.  The  experimental  observations  reported 
here  increase  our  knowledge  of  this  behavior  and  the 
physical  model  both  increases  our  understanding  of  the 
phenomenon  and  provides  us  with  a reliable  method  of 
predicting  the  compressibility  of  wet  snow  under  a 
variety  of  conditions.  The  model  is  in  fair  agreement 
with  the  experimental  observations  over  a wide  range 
of  conditions. 
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